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Abstract. There have been exciting advances recently in finite temperature 
calculations with three quark flavours and with non-zero chemical potential. The 
role of improved actions is explained and recent results from the Bielefeld group and 
MILC collaboration with improved Kogut-Susskind quarks are presented. Three new 
approaches to finite chemical potential are discussed. 
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1. Introduction 

I have been asked by the organizers to review the status of lattice QCD calculations at 
non-zero temperature with three quark flavours. The exciting news here is that due to 
improvements in algorithms, we are now able to carry out calculations with light up and 
down quarks considerably lighter than the strange quark. Before presenting results, I 
give a brief introduction to lattice calculations, and then I will outline how algorithms 
have been improved. After presenting numerical results, I will discuss briefly the recent 
significant advances in dealing with finite chemical potential. (I ran out of time during 
my talk and did not discuss this at the conference.) 

2. Introduction to Lattice Calculations 



To carry out a lattice simulation we must select certain parameters, namely the lattice 
spacing (a) or gauge coupling (J3), a fixed grid size (N% x N t ) where N s and N t are 
the space and time dimensions of the grid, respectively, and quark masses (m u ^, m s ). 
There are also certain parameters related to the algorithm. The parameters detailed 
have physical meaning, and each choice can lead to a systematic error that must be 
controlled. 

To deal with the non-zero lattice spacing, we must take the continuum limit. 
Similarly, we must take the infinite volume limit to deal with the finite fixed grid. 
Finally, we must extrapolate to light quark mass for the up and down quarks. This is 
a practical issue as it is too expensive to do the calculations with the physical up and 
down masses. On the other hand, we can work at the physical s quark mass. For finite 
temperature calculations, the temperature is given by T = l/(N t a), so we use grids with 
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N t < N s . Typically, N s is two or three times N t to provide a reasonably large physical 
volume. 

Also of physical relevance (especially for the RHIC program) is nonzero chemical 
potential /i. Until recently, it has been nearly impossible to carry out calculations with 



3. How have algorithms been improved? 

The first numerical lattice calculations used the Wilson (plaquette) gauge action and 
either Wilson or Kogut-Susskind (KS or staggered) quark actions. The Wilson quark 
action is designed to solve the fermion "doubling" problem. So is the KS approach; 
however, it maintains enough degrees of freedom for four quark fields. There remains a 
U(l) chiral symmetry and a single Goldstone pion. The rest of the flavour symmetry 
is broken for non-zero lattice spacing. Thus, the other pion states are heavier than the 
Goldstone state because of flavour (or what is now frequently called taste) symmetry 
breaking. 

Starting in the mid-1980s, the Symanzik improvement program P began to be 
applied to calculations to improve the scaling properties of the theory. For Wilson 
type quarks the improvement called clover quarks was introduced by Sheikholeslami 
and Wohlert [2|. Naik introduced a new three link term for KS quarks in 1989 j^j. 
An improved gauge action was applied in 1995 j3j. Soon thereafter calculations were 
done with the Naik term and an improved gauge action [S] and the importance of 
"fattening" the action to reduce taste symmetry breaking was tested and understood 
[6-10]. About the same time, the P4 action was introduced because of it's improved 
rotational symmetry 

All of these improvement schemes gous to higher order methods in 

numerical analysis, but with the twist of applying to QFT. They exact a computational 
toll from having a more complicated action, but they more than make up for that 
by reducing systematic errors and allowing a larger lattice spacing to be used. Since 
computational requirement goes like a large power (7-9) of the inverse lattice spacing, 
it really reduces the cost to use a larger lattice spacing. 

More recently, domain wall fermions and the overlap method were developed to 
better control chiral symmetry. They are not yet extensively used for thermodynamics, 
but see reference jTSj. 

We shall concentrate on the improvement program for staggered quarks. The 
Bielefeld group is the center of activity for the P4 action that maintains rotational 
invariance of the free quark propagator to order p , 
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The MILC collaboration uses the "Asqtad" action jHJElEn]- Their gauge action 
includes plaquette, 1x2 and bent 6-link terms. Their quark action includes the links 
shown below and the 3-link Naik term. The three, five and seven link terms are known 
as fat links and are necessary to reduce taste symmetry breaking. The full set of fattened 
links is pictured in figure The weights of each term are not shown. Each diagram 
represents a term in the quark action of the form x( x )V(x,x + fi)x( x + /*)> where V 
is the product of links along the path. With coefficients of the various terms set using 
tadpole improvement, the errors of this action are of order a 2 a and a 4 . 




Figure 1. Terms used to suppress flavour symmetry breaking. The final five link 
path was introduced by Lepage to correct the small momentum form factor. 

It is easy to see why improved actions should have better high temperature behavior 
than the Wilson or KS actions. In figure we can examine the energy and pressure 
for free massless quarks on a lattice with N t time slices [H]. The continuum limit 
corresponds to Nt —>■ oo. One immediately sees that both Naik and P4 actions much 
more rapidly approach the continuum limit than either Wilson or KS quarks. At N t = 4, 
Naik is better for pressure and P4 is better for energy density. For N t — 6, the situation 
is reversed. For N t > 8, both improved actions are quite close to the continuum limit, 
with P4 a little closer. 

4. Some recent physics results 

We will now concentrate on results with 3 degenerate or 2+1 flavours of quark. Some 
of the issues to be addressed are: a) what is the phase diagram of QCD?, b) what is 
the transition temperature?, c) where does the the physical point in the phase diagram 
lie?, d) what is the equation of state for QCD? 
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Figure 2. The energy density (left) and pressure (right) of free massless fermions as 
function of temporal lattice size N t . 
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Figure 3. From F. Karsch [T5] , 



To find the phase diagram, we must find the transition temperature for a variety 
of quark masses, so it would be logical to address (b) before (a); however, presenting 
the phase diagram first helps to orient us, so in figure El I present a nice summary plot 
prepared by Karsch of the QCD phase diagram in the m u ^-m s plane [To] . 

Symmetry considerations help us to identify the order of the transition on the 
corners of the diagram [Ej. In the upper right hand corner we have pure gauge theory 
with a deconfinement transition at T ~ 270 MeV. In the lower- left corner of the diagram, 
we expect a first order transition for three massless degenerate flavours. Chiral symmetry 
is broken at low temperature and restored at high temperature. We expect that as 
the quark masses are increased, the transition will weaken as we approach a second 
order line. When the strange quark is very heavy, but u and d are light, we expect a 
second order transition for massless quarks. Karsch's estimates of the chiral transition 
temperatures in the corners are shown in the diagram. 

Some recent results from the MILC collaboration with the Aqstad action and 
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Figure 4. ("ipip) vs. temperature for various light quark masses. Bursts are 
extrapolations to zero quark mass. These graphs are updated from the ones in 
reference |17j . 
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Figure 5. (ipip) for three degenerate quarks with Nt =8. 



either 2+1 or 3 degenerate flavours indicate that above 175-200 MeV (ipip) extrapolates 
to as m q — > 0. This is evidence that chiral symmetry is restored at higher temperature. 
For 2+1 flavours, with N t = 6, the chiral symmetry restoration starts at about 200 MeV. 
The runs with N t = 8 continue as the lightest quark mass runs are not completed. Here 
it appears the transition may occur at about 175 MeV. The difference can be a finite 
lattice spacing effect. Calculations have also been done with three degenerate quarks. 
The lightest quark mass explored by MILC is 0.4 m s where m s is the strange quark 
mass. Each of the individual curves looks quite smooth, so it is not clear that any 
of these masses are within the first order transition region; however, if the curves are 
extrapolated to zero quark mass, there appears to be chiral symmetry restoration around 
185 MeV. 
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Figure 6. Ratio of transition temperature to vector (rho) mass vs. squared ratio of 
pseudoscalar to vector masses |18) . Results are shown for various numbers of dynamical 
quarks. For Nf = 2 results with an unimproved action are shown for comparison. 

Turning now to the P4 action, Karsch, Laermann and Peikert have done an 
extensive study with N t = 4 A graph summarizing their results for the ratio 

of T c to the vector meson mass is show in figure H3 For Nf = 2 they show both the 
P4 action and the standard KS action results. Clearly, the improved action gives a 
higher value for T c . With three degenerate flavours, the transition temperature is lower 
than for Nf = 2 (considering only improved action results). One combination of quark 
masses is available for Nf = 2 + 1; however, the result looks very much like the Nf = 2 
result. This could be because the strange quark mass is so much larger than its physical 
value. The Bielefelders estimate that in the chiral limit T c = 173 ± 8 MeV for Nf = 2 
and 154 ± 8 MeV for N f = 3. 

We see that the two groups do not seem to be in close agreement on the transition 
temperature for Nf = 3. A careful comparison would be useful. It is possible that the 
difference is a finite lattice spacing effect. 

We would like to find the line in the phase diagram separating the first order 
transition region for very light quarks from the crossover region for heavier quarks. 

MILC is exploring a horizontal line at fixed m s and the diagonal with m s = m u d- 
The results shown for {^pip) are rather smooth curves except perhaps for the lightest 
quark mass studied on the 16 3 x 8 lattice. MILC does not claim to have worked at light 
enough quark mass to be in the first order transition region. Two other groups have 
also studied this question. Karsch, Laermann & Schmidt pH] and Christ & Liao [20] 
have used light quarks with N t = 4 and the unimproved action. Both groups are now 
using the Binder cumulant method to determine the order of the transition. In figure 
the Binder cumulant is shown as a function of quark mass for three different spatial 
volumes. The lines should cross at the critical point, which appears to be about 0.035 
in lattice units. 




Figure 7. Binder cumulant from Karsch, Laermann & Schmidt [T2|. Results from 
different volumes should cross at the critical point separating the first order region 
from the smooth crossover. 



The two groups estimate the pseudoscalar mass at the critical point to be 270 and 
290 MeV. Karsch et al. have some results with the P4 action for which M ps « 192(25) 
MeV. Clearly, given the difference between the unimproved and P4 actions, more work 
to control the cutoff effects is needed to capitalize on these exploratory works. Also, as 
interesting as is the three degenerate flavour case, we are most interested in m s > m U( i. 

The pressure of the quark-gluon plasma is of great interest. Using the P4 action, 
Karsch, Laermann and Peikert [21] have done an extensive study for N t = 4. Good 
control of flavour symmetry is important near T c because there is a system of light 
hadrons, and flavour symmetry breaking distorts the spectrum of the light hadrons. 
Some fattening is used in this calculation to improve flavour symmetry, but it would 
be necessary to repeat the calculation with large N t , i.e., smaller lattice spacing to 
get reliable results near T c . However, far above T c , these results may be indicative 
of the continuum limit. Figure |H1 shows the pressure using the P4 action for various 
combinations of quark mass and flavours. Indicated by arrows are the free quark values 
at very high temperature. In the second graph, each pressure curve is normalized by 
the corresponding value. A comparison with unimproved action results and an estimate 
of the continuum limit at high temperature is in figure 01 



5. Strange Quark Content of QGP 

MILC has looked at various quark number susceptibilites They have been related 
to event-by-event fluctuations in heavy ion collisions in work by B. Muller [22] using the 
fluctuation-dissipation theorem. 

(tQ 2 )«y^t l = x Q (T,» Q = 0), (2) 
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Figure 9. Pressure estimate by Karsch, Laermann and Peikert [21] of continuum limit 
for T > 2T C . 



We define: 

T d 2 log Z 



Xij 



(3) 



V s dfjLid/jLj 

Xsing ^Xuu ^Xud (4) 

and triplet quark number susceptibility 

Xtrip = 2Xuu — ^Xud ■ (5) 

As seen in the second part of figure ITUl the difference between Xsing an d Xtrip is most 
pronounced around the transition temperature. MILC also has results for the strange 
quark susceptibility. It would be useful to try to turn these results into a statement 
about the strangeness content of the QGP in a way that is testable at RHIC. 

6. Nonzero Chemical Potential 



Due to lack of time, this material was not presented in my talk. That was unfortunate, 
because very significant progress has been made in algorithms, and we can expect to 
see concomitant advances in calculations soon. 
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Figure 10. Left: The triplet and strange quark number susceptibilities for Nf = 2 + 1 
with m u d — 0.2 m s on 12 3 x 6 and 16 3 x 8 lattices Right: The difference 

between singlet and triplet quark number susceptibility for Nf = 3 with quarks of 
mass niq = 0.4 m s , on 8 3 x 4, 12 3 x 6 and 16 3 x 8 lattices. 



Calculations with finite chemical potential are important because experiments at 
heavy-ion colliders such a RHIC start with nuclei, not anti-nuclei so there is a quark 
chemical potential \x of about 15 MeV |2S]- However, finite chemical potential is difficult 
because the action is no longer real so it is no longer possible to implement importance 
sampling in the usual way. Recent progress has been made using three new methods. 

Fodor and Katz introduced a multiparameter reweighting technique [21]. The key 
here is that the gauge coupling is shifted from (3 to /3q in order to maximize the overlap 
of the ensemble generated with \i = with the desired \i ^ ensemble. This is detailed 
in (0) and illustrated in figure [TT] [23] . 



Z(J3,fi,m) = J T)Uexp[-Sg(J3)]detM(n,m) 

= J VUexp[-S g ({3 )]detM(ij = 0,m) (6) 

A second method introduced by a Swansea-Bielefeld group uses a Taylor series 
expansion in fi/T of both the reweighting factor and observable variables |26| . Thus, 
reweighting is avoided but the chemical potential must be small as only the leading 
order term is available. 

de Forcrand and Philipsen use analytic continuation from imaginary \i |28| l2"Tj . 
When the chemical potential is imaginary, the action is real and importance sampling 
works. The challenge is to do the analytic continuation. It will be interesting to see 
which of these approaches is best. We are certainly learning more about QCD with 
/i 7^ than we had previously been able to. 
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Figure 11. Best weight lines in fj, — j3 plane |25| . 
7. Prospects 

I hope I have convinced you that there are many interesting calculations being done 
at finite temperature with Nj = 2 + 1 and 3. However, a great deal remains to be 
done: we need better control of systematic errors, including the continuum limit and 
wider coverage of the phase diagram. I have concentrated on results with improved KS 
type quarks, but these are not the only methods for quarks. For example, improved 
actions for Clover quarks are being pursued by the CP-PACS/JLQCD collaborations 
[2TI] . Dynamical quark calculations with domain wall or overlap quarks and calculations 
with chemical potential are just in their infancy. Although we lattice practitioners 
have learned a great deal recently, there is much we can and shall do to improve our 
calculations. 
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